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ABSTRACT 


A new analytical model for simulating chemoviscosity of thermosetting 
resin has been formulated. The model is developed by modifying the 
well-established Wi 11 lams-Landel -Ferry (WLF) theory in polymer rheology for 
thermoplastic materials. By assuming a linear relationship between the glass 
transition temperature Tg(t) and the degree of cure a(t) of the resin system 
under cure, the WLF theory can be modified to account for the factor of 
reaction time. Temperature-dependent functions of the modified WLF theory 
constants C^T) and C 2 (T) were determined from the isothermal cure data of 
Lee, Loos, and Springer for the Hercules 3501-6 resin system. Theoretical 
predictions of the model for the resin under dynamic heating cure cycles were 
shown to compare favorably with the experimental data reported by Carpenter. 
This work represents progress toward establishing a chemoviscosity model which 
is capable of not only describing viscosity profiles accurately under various 
cure cycles, but also correlating viscosity data to the changes of physical 
properties associated with the structural transformations of the thermosetting 
resin systems during cure. 
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CHEMOVISCOSITY MODELING FOR THERMOSETTING RESINS - II 


I. INTRODUCTION 

There are two different basic approaches in chemovi scosity modeling of 
thermosetting resin. One approach, typically represented by the work of 
Roller [1], is to empirically formulate a rate equation which relates the 
change of chemovi scosi ty with reaction time. For a resin system cured under a 
dynamic heating cure cycle (as commonly encountered in autoclave processing 
for composite materials), such an approach usually yields a model which is 
inadequate in describing accurately the nonlinearity of chemovi scosi ty as a 
function of reaction time. As Tajima and Crozier [2] pointed out, such a 
modeling approach renders itself to the limitation of batch-specific. The 
model parameters cannot be readily related to the chemical and rheological 
properties of the reacting system as well. 

The second approach is based upon a modification to the well-established 
vi scosi ty-temperature relationship existing in polymer rheology for 
thermoplastic materials. The parameters in such an equation can be expressed 
in terms of polymerization kinetics, and the chemovi scosi ty profiles as a 
function of reaction time can then be modeled for thermosetting resin systems. 

The applicability of the modified Wi 1 1 i ams-Landel -Ferry (WLF) theory [3] 
in chemovi scosi ty modeling for thermosetting resin has been studied by Tajima 
and Crozier [2,4], Apicella et al . [5], and Hou [6] among others. It has been 
extensively documented in the literature that temperature-dependent viscosity 
of thermoplastic material can be accurately described by the WLF equations 



within 100°K above its glass transition temperature (T g ). For a 
thermosetting resin system under cure, it is reasonable to assume that Tg(t) 
of the material is always within 100°K lower than the cure temperature. The 
WLF equation should then be applicable to every instant state of the resin 
system under curing. The central theme in modifying the WLF theory to 
describe the chemoviscosity of thermosetting resin is to take into account the 
reaction time factor. 

In an earlier paper [6], two assumptions were proposed, namely, that the 
rate constant of reaction, kj, at any temperature T, is (i) diffusion 
controlled and is, therefore, inversely proportional to the viscosity 
n-j-(t) of the reactive medium, and ( 11 ) directly proportional to the rate 
of change of glass transition temperature T g (t). The modified WLF equation 
becomes a first order ordinary nonlinear differential equation. Numerical 
solutions have also been shown to compare favorably with the experimental 
results for several thermosetting systems under isothermal and dynamic heating 
cure conditions. It has been concluded that the flexibility demonstrated by 
the modified WLF theory [6] can be conveniently exploited to establish an 
analytical model with high degree of accuracy for the chemoviscosity of any 
thermosetting resin system under various cure cycles. The physical 
significances of the material parameters selected for the model were, however, 
difficult to extract for the particular resin system under investigation. 

As a next step toward the improvement of the chemoviscosity modeling for 
a thermosetting resin system, a new modification of the WLF theory is proposed 
in this report. Comparisons between the theoretical predictions and 
experimental data for chemoviscosity of a resin system (Hercules 3501-6) under 
both isothermal and dynamic heating cure conditions will be discussed. 
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List of Symbols 


^ 1 ’ A 2 *^3 
3 9 b , c 9 d 
B 

C ! ( T ) , C 2 (T) 


DSC 


AE 1 , AE 2 , aE 3 


K i ,K 2 - 


R 


T 

T C 

Tg(t) 


Tg(°°) 

T° 

g 

t 

tc 

tf 

At 

WLF 

Tlj(t ) 

cc(t) 


Material parameters (min -1 ) specified in Table 1 
Constants used in Eq. (5) 

Material parameter specified in Table 1 
Temperature-dependent WLF constants defined by Eqs. (7) and 
( 8 ) 

Differential Scanning Calorimeter 
Activation energies (J/mol) specified in Table 1 
Material constants (min -1 ) defined in Eqs. (3) and (4) 
Universal gas constant (=1.987 Cal /gm-mole/°K) 

Temperature (°K) 

Cure temperature (°K), T c = Tg(<») 

Glass transition temperature of the resin system 
under cure (°K) 

Glass transition temperature of the resins system at 
completion of cure (°K) 

Glass transition temperature of the uncured resin g 
system (°K) 

Cure time (minutes) 

Time (minutes) defined in Eq. (6) 

Time required for completion of cure (minutes) 

Time interval (minutes) used in solving Eq. (9) 

Willi ams-Landel -Ferry 

Chemoviscosity at temperature T (poises) 

Degree of cure of the thermosetting resin during cure 
as determined by DSC thermal analysis 
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II. THEORY 


The Wi Hi am-Landel -Ferry (WLF) theory [3] is given by 



) - 


-Ci(T-T q ) 

C 2+ T-T g 


( 1 ) 


for Tg < T < T g + 100° K. 

C x and C 2 are two material constants, and n-p represents viscosity at 
temperature T of the given thermoplastic material which possesses a glass 
transition temperature Tg. The normal use of the WLF equation for 
thermoplastic materials requires that Tg be constant while the temperature T 
is varied for the specific polymer under study. However, during cure of 
thermosetting resins, the monomers are initially polymerized and crosslinks 
are formed later. This is a system where Tg(t) is changing and the curing 
temperature T is held constant (in an isothermal cure case, for example). The 
glass transition temperature T g rises continuously and may eventually 
approach the curing temperature. Over the entire curing cycle, the material 
structure actually undergoes continuous phase transformations from a low 
molecular weight liquid to a high molecular weight polymeric melt, and 
eventually transforms to a crossl inked network. It is reasonable to assume 
that Tg ( t ) of the material is always lower than the cure temperature T, and 
that (T-Tg(t)) is always within 100°K. The WLF theory (Eq. (1)) should then 
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be applicable to ail of the different polymeric structure phases during cure. 
Before Eq. (1) can be applied to describe the chemovi scosity of the curing 
resin, however, modifications need to be made in order that reaction time, t, 
of the thermosetting resin system is properly taken into account. This can be 
accomplished as follows: 

A simple assumption that Tg(t) is a linear function of the degree of 
cure a(t ) of the thermosetting resin can be made, which gives: 

yt> - T g° * <y T g> «<‘> < 2 > 

where T° = ( 0) is the glass transition temperature of the unreacted resin 

system, T c is the cure temperature in °K, and t is the reaction time in 
minutes. Equation (2) should satisfy the boundary conditions, namely, 

T = T°, a = 0 at t = 0 

g g 


and 


Tg = T c , a = 1 at t = t f 

where tf is the time required for completion of cure. The degree of cure 
a(t ) for the resin system can be determined by means of thermal analysis. It 
is assumed that the degree of cure at time t is equal to the fraction of heat 
released as measured by a DSC up to time t for the thermosetting material 
under cure. The total heat of reaction can be determined for the experiments 
using dynamic temperature scanning. The degree of cure a(t) can then be 
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determined for different isothermal cure conditions [7,8], The a's thus 
determined are functions of curing temperature and time. 

Equations (1) and (2) constitute a modified WLF theory for simulating 
chemoviscosity of thermosetting resin. Among the five material parameters, 
and C 2 can be determined as a function of temperature T by fitting with 
isothermal cure n(a) data, n T = 10 13 poises is given by Gordon and 

'g 

Simpson [9] for many polymeric materials at T = T q , and T° can be determined 

y 9 

by DSC thermal analysis. In an isothermal cure case, T c of Eq. (2) is equal 
to cure temperature. In the dynamic heating case, however, T c is difficult 
to specify exactly and is the only adjustable variable in the present modified 
WLF theory. 
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III. EXPERIMENTAL DATA 


For the purpose of evaluating model predictions of this modified WLF 
theory, experimental data on Hercules 3501-6 epoxy resin reported in the 
literature by Lee, Loos, and Springer [8] and Carpenter [10] were selected. 
These data will be summarized briefly. 

Degree of Cure a(t,T) 

Lee et al . [8] performed thermal analyses on the epoxy resin system using 
a Differential Scanning Calorimeter (DSC). The degree of cure a and the rate 
of degree of cure da/dt were determined from the results of the isothermal 
scanning experiments. Eight different temperatures were selected, and the 
results are reproduced in Figure 1. It is noted from the figure that both 
higher values of the degree of cure and rate of degree of cure can be obtained 
at higher cure temperatures. The rate of degree of cure is a decreasing 
function of cure time t, and eventually levels off at a low rate value (< 0.1) 
when a high degree of cure is reached. The following equations were found to 
describe these temperature-dependent data accurately: 


= (K^Kjct) (1-cO (B- a ) for ct^O.3 


(3) 


da 

St = 


(l-a) 


for a > 0.3 


(4) 


where 
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K l = ^1 ex P (- AE i/RT) 
K 2 = A 2 exp (-aE 2 /RT) 
K 3 = A 3 exp (-aE 3 /RT) 


Values of constants in Eqs. (3) and (4) are shown in the Table below. 
Comparisons between the measured values of da/dt and the values calculated by 
Eqs. (3) and (4) are also shown in Figure 1. 

Table 1. The values of the constants of Eqs. (3) and (4) 
for Hercules 3501-6 Resin System [8] 

B = 0.47 

A x = 2.101 x 10 9 min* 1 
A 2 = -2.014 x 10 9 min" 1 
A 3 = 1.960 x 10 s mm- 1 
AE! = 8.07 x 10 4 J/mol 
aE 2 = 7.78 x 10 4 J/mol 
aE 3 = 5.66 x 10 4 J/mol 


Eqs. (3) and (4) can be integrated provided that the temperature is kept 
constant. The following relationships between a and t are obtained: 

a K 2 

t = Ln (1 + y- a) - b Ln (1-a) - c Ln (1- f-) (5) 


t 


K, Ln ( 0.7* + t c 


( 6 ) 
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where 


a = K 2 (B-l)/d 
2 

b = - ( K !+K 2 ) /d 
c = (K 1 +K 2 B)/d 

d = (K^+K 2 ) [K jB 2 +( K j -K 2 ) B-K x ] 
and t c is the time given by Equation (5) for a = 0.3. 

Chernov i scosi ty n(t,T) 

Chemovi scosi ty of the resin under an isothermal cure condition were 
measured by Lee et. al . [8] on a Rheometncs parallel disc and plate 
mechanical spectrometer Model 605. The radius of the disc was 25 mm and the 
gap between the disc and the plate was 0.5 mm. A frequency of 1.6 Hz was 
employed. For a given cure temperature, the n(t) measured can be related to 
a(t) obtained from DSC thermal analysis. The results are plotted in Figure 2 
for three different isothermal cure cases. It is again noted that at a higher 
cure temperature, higher values of the degree of cure can be reached. For a 
given a, the decrease in viscosity is attributed to the temperature effects. 

Chemovi scosi ties of the 3501-6 epoxy resin cured under dynamic heating 
conditions were reported by Carpenter [10]. Carpenter measured chemovisco- 
sities at different rate of increase in temperatures and reported the measured 
values as a function of cure time. The experimental data shown here are for 
resins cured at constant heating rates of l°K/min (figure 5), 2.17 and 5°K/min 
(Figure 6) and a cure cycle (Figure 8) similar to the one typically found in 
autoclave processes for composite materials. 


9 



The procedure used for model evaluations will be as follows: 

1. Model parameters of Eqs. (1) and (2) will be determined first by the 
experimental data of Lee et al . [8] under isothermal cure conditions. 
Temperature-dependent material constants C^T) and C 2 (T) can thus be 
determined. 

2. The chemoviscosity model n(t,T) established by isothermal cure data under 
procedure 1 will then be made to follow various dynamic heating cure 
cycles T{t). The viscosity values thus calculated will be compared with 
the experimental data reported by Carpenter [10]. 
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IV. MODEL PREDICTIONS AND DISCUSSION 


/ 


Isothermal case 


Gordon and Simpson [9] pointed out that glass-forming substances 

generally have the same shear viscosity of n T 3 1Ql3 poises at their 

9 

Tg, and this has also been confirmed for polymer systems. By specifying T° 

T° = 300°K (27°C) in Eq. (2), the model (i.e., Eqs. (1) and (2)) can then be 
fitted to the experimental data n(a) as shown in Figure 2 provided that proper 
values of material constants C : and C 2 are selected. C x and C 2 determined for 
different temperatures are tabulated below, and the model predictions are 
drawn as solid lines in the figure. 


Table 2. Values of material constants C L and C 2 
in the modified WLF Theory for Hercules 
3501-6 resin system as determined by the 
isothermal cure data. 


T C («K) 

Ci 

c 2 (°k) 

375 

21.1 

53.0 

400 

17.1 

33.6 

425 

18.5 

46.0 


The model is capable of describing the nonlinear relationship between 
viscosity, n» and the degree of cure, a, as shown in Figure 2, especially for 
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the case of higher cure temperature. This is considered a significant 
improvement over the empirical linear model used by Lee, Loos, and Springer 
[8]. An arrhemus type of plot is arbitrarily chosen for Cj and C 2 as shown 
in Figure 3. The data is rather scattered. Due to the (unfortunate) 
limitation on size of the experimental data and the wide range in temperatures 
covered by the dynamic heating cure experiments (as will be discussed later), 
a linear least squares fit was chosen. The straight lines are represented, 
respectively, by: 

log C x = 186.02 (1/T) + 0.8084 (7) 

and 

log C 2 = 212.16 (1/T) + 1.106 (8) 

It is noted that in the analysis described above, the glass transition 
temperatures Tg(t) for the varying molecular structures of the resin cured 
at different tempreatures are chosen as reference temperatures. C 1 and C 2 are 
therefore treated as "universal constants" in the present context. Values of 
C : and C 2 reported in Table 2 are close to the values of 17.44 and 51.6, 
respectively, once reported by Williams et al . [3] as the universal constants 
for all polymeric materials sytems. Eqs. (7) and (8) will still be used, 
however, to account for temperature-dependences of C x and C 2 in chemovi scosity 
simulation under dynamic heating cure conditions. 

The theoretical assumptions as described above, which yield close values 
for the WLF constants and C 2 , respectively, among three different 
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isothermal curing temperatures, are also supported by the parallelism among 
experimental data observed in Figure 2 for a < 0.3. Data of higher values of 
both n and a are difficult to be measured, and could be considered less 
accurate. 

Dynamic Heating Case 

In the case of dynamic heating cure conditions, integrations of Eqs. (3) 
and (4) are very difficult to perform. If one assumes that within any time 
interval At = t 2 - t x , the temperature T c is uniform, and is equal to 
( T i +T 2 ) / 2 , where Tj and T 2 are temperatures at t l and t 2 respectively as 
defined by the cure cycle, Eqs. (3) and (4) become: 

At = 1 2 ~ 1 1 

= f U 2 .T c ) - f(c4,T c ) (9) 


where 


f 


(a,T c ) = ^ Ln (1 



a) - b Ln ( 1 -a) - C Ln 


(1 


- f-)i for a _< 0.3 


and 


f (a * T c ) = _ r; Ln( V^ ) + t c 


for a > 0.3 
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The parameters a, b, c, and t c are the same as defined before in Eqs. (5) 
and (6). The value of the increased degree of cure a 2 in Eq. (9) can be 
calculated by the Generalized Newton's Method. Details of the solution 
schemes are summarized in the Appendix. 

Values of ct(t) determined by the scheme described above depend heavily on 
the time interval At chosen. Theoretically more accurate a(t) can be 
determined numerically when a smaller At is used. On the other hand, a 
smaller At could increase dramatically the numerical round-off errors. A 
series of At values were selected, and a(t) was calculated from Eq. (9) for a 
cure cycle with a constant heating rate of 5°K/min. The results are plotted 
in Figure 4. It can be seen that numerical round-off errors become dominant 
when At < 0.05 min. It is also noted that numerical values of a(t) calculated 
for At = 0.5 and 1.0 minute, respectively, are almost indistinguishable. 
Consequently, At = 0.5 minute was chosen for all calculations discussed below. 

Chemovi scosities of the Hercules 3501-6 epoxy resin system under dynamic 
heating cure conditions were reported by Carpenter [10], Data for these 
constant rates of heating experiments are reproduced in Figures 5 and 6. The 
heating rates are 1.0, 2.17, and 5 °K/min respectively. All cure cycles 
started from 300°K (27°C). Temperature-dependent Cj(T) and C 2 (T) as defined 
by Eqs. (7) and (8) are followed in the simulation. Values of the rest of the 
parameters used in the model are tabulated below: 
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Table 3. Values of model parameters used in simulating 
chemov i scosity under rate of heating cure 
conditions 


Heating rate 
(° K/min) 

ny (poises) 

g 

T° (° K) 

T c (°K) 

1 

10 13 

300 

415 

2.17 

10 13 

300 

430 

5 

10 13 

300 

450 


The model predictions are shown as solid curves in Figures 5 and 6. It is 
noted that the theoretical calculations compared favorably with the 
experimental data for all three cases. 

Theoretically, T c is equal to Tg( = ) at T + « or a > 1 (see Eq. (2)). 
Under the dynamic heating cure conditions, however, T c is difficult to 
specify. Values of T c (in Table 3) in the above calculations were 
selected arbitrarily such that model predictions are in reasonably good 
agreement with the experimental data. Nevertheless, higher values of T c 
(Table 3) resulting from the simulation of higher rates of heating 
experiments, are in agreement with physical observations for thermosetting 
resins . 

The calculated degree of cure a(t) for all three cases is plotted in 
Figure 7. These calculated a(t) are based on the experimental data for the 
resin under isothermal cure conditions [8 ] , and are independent of the 
arbitrarily selected values of T c discussed above. Values of da/dt go 
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through maxima for all three cases. The arrow marks shown in the figure 
denote the times when chemoviscosities approach infinity during cure. These 
arrow marks are located approximately at a(t) = 0.75 _+ 0.05, and are near the 
inflection points of the curves as well. Physical significances of these 
observations are unclear at the present time, and further research in the 
future is required. 

Chernov iscosity was also measured by Carpenter [10] for the 3501-6 resin 
system under a typical cure cycle used in autoclave processing for composite 
materials. The experimental data is reproduced in Figure 8. Values of model 
parameters used are the same as those included in Table 3 and Eqs. (7) and 
(8), except that T c = 375° K. The theoretical predictions, shown as solid 
curves, are seen in reasonably good agreement with the experimental data. The 
model, however, fails to describe the chemoviscosity "shoot-up" near the end 
of the cure cycle. The reasons are unclear at the present time. It could be 
due to the inadequate definitions of C X (T) and C 2 (T) by Eqs. (7) and (8); or 
due to different reaction mechanisms appeared in the high curing temperature 
regime. Calculated values of a(t) are plotted in Figure 9. The arrow mark 
shown denotes the time when chemoviscosity "shoot-up" occurs. Features of the 
a(t) curve discussed above for isothermal cure cases are again noticeable in 
this figure. 

Temperature ranges covered by Carpenter [10] in dynamic heating cure 
experiments, as shown by dotted lines in Figure 3, far exceed the ranges 
covered by Lee et al . [8] in isothermal cure experiments (shown by solid lines 
in the figure). Accurate definitions of the crucial temperature-dependences 
of C^T) and C 2 (T) needed for the model thus suffer from serious data 
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extrapolation errors. The results of the present investigations, although 
encouraging, can only demonstrate the feasibility of the approach of modifying 
the WLF theory in chemoviscosity modeling for thermosetting resin. Detailed 
studies such as the precise definitions of Tg(a) (Eq. (2)), and the 
functions of C^T) and C 2 (T) over the full temperature range by experimental 
measurements for the resin system under investigation have to be pursued in 
the future. 
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V. CONCLUSIONS 


A new analytical model, based upon modifications of the WLF theory, has 
been formulated. The model assumes that the glass transition temperature 
Tg(t) of the resin system under cure is linearly related to the degree of 
cure a(t) determined by DSC thermal analysis. Theoretical predictions of the 
chemov iscosity compare favorably with the experimental data reported in the 
literature for both isothermal and dynamic heating cure conditions. 
Considering that data are available only in the relatively narrow temperature 
range used in isothermal curing experiments by Lee, Loos, and Springer [8], 
the agreement between theoretical predictions and experimental data reported 
by Carpenter [10] for dynamic heating cure cycles are particularly 
encouraging. Due to the limitations in size of the data available, the 
successful analysis reported here has established the feasibility of the 
proposed approach which simulates chemoviscosity for thermosetting resin by 
modifying the widely used WLF theory in polymer rheology. 
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Appendix - Generalized Newton's Method in Solving <* 2 from Eq. (9) 


Eq. (9) can be rearranged to become: 

g( ct 2 ) = f (« 2 » T c ) - f (oti,T c ) - At c = 0 


(a) 


and 


g’(a 2 ) = (b) 

where g'(a 2 ) = dg/da 2 , and f'(a 2 ,T c ) = df/da 2 . The Generalized Newton's 
Method sets up the following recursion formula to solve for a 2 in Eq. (a): 

a 2 ^=a 2 -w.g(a 2 )/g'(a 2 ),n=0,l,2... (c) 

provided that g'( a 2 )/*0; where w is called the over-relaxation factor. A good 
choice for w can be determined, in general, only by experimentation. A choice 
of w different from unity often can increase the convergence rate of the 
routine appreciably. In this study, w = 0.1 was used for all calculations. 
Iterations are continued until |a n+ ^ - a n |< 10 -6 . 
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Rate of degree of cure vs. degree of cure reported by Lee, Loos, 
and Springer [8] for Hercules 3501-6 resin system under isothermal 
cure cycles. Solid lines were calculated using Eqs. (3) and (4). 

Chernov iscosity as a function of degree of cure for Hercules 3501-6 
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experimental data reported by Lee, Loos, and Springer [8], Solid 
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by the isothermal cure data of Lee et al . [8] for the Hercules 
3501-6 resin system. 

Simulated degree of cure data for the Hercules 3501-6 resin system 
cured under constant rate of heating condition, 5°K/min. Various 
time intervals At were used in solving Eq. (9). 

Chemoviscosity of the Hercules 3501-6 resin system cured under 
constant rate of heating condition of l°K/min. Circles denote 
experimental data by Carpenter [9]. Solid curve denotes the model 
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Chemoviscosity of the Hercules 3501-6 resin system cured under 
constant rates of heating conditions of 2.17 and 5°K/min, 
respectively. Circles denote the experimental data by Carpenter 
[9]. Solid curves denote the model predictions of Eqs. (1) and 
( 2 ). 

Simulated degree of cure for the Hercules 3501-6 resin system cured 
under constant rate of heating conditions indicated. 

Chemoviscosity of the Hercules 3501-6 resin system cured under the 
cure cycle shown. Circles denote the experimental data by 
Carpenter [9]. Solid curve was calculated by the model, Eqs. (1) 
and (2). 

Simulated degree of cure data for the Hercules 3501-6 resin system 
cured under the cure cycle included. 
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RATE OF DEGREE OF CURE, da/dt 



DEGREE OF CURE, a 


Figure 1. Rate of degree of cure vs. degree of cure 

reported by Lee, Loos, and Springer [8] for 
Hercules 3501-6 resin system under 
isothermal cure cycles. Solid lines were 
calculated using Eqs. (3) and (4). 
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Figure 2. Chemovi scosi ty as a function of 

degree of cure for Hercules 3501-6 
resin system under isothermal cure 
cycles. Symbols represent 
experimental data reported by Lee, 
Loos, and Springer [8]. Solid 
lines were calculated by the model, 
Eqs. (1) and (2). 


24 



in, °k x w 


Figure 3. Temperature-dependences of the WLF 

constants and C 2 determined by the 
isothermal cure data of Lee et al . [8] 
for the Hercules 3501-6 resin system. 
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Figure 4. Simulated degree of cure data for the Hercules 
3501-6 resin system cured under constant rate 
of heating condition, 5°K/min. Various time 
intervals At were used in solving Eq. (9). 





Figure 5. Chernov i scosi ty of the Hercules 3501-6 resin 
system cured under constant rate of heating 
condition of l°K/min. Circles denote 
experimental data by Carpenter [10]. Solid 
curve denotes the model predictions of Eqs. 
(1) and (2). 
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Figure 6. Chemovi scosity of the Hercules 3501-6 resin 
system cured under constant rates of heating 
conditions of 2.17 and 5°K/min, respectively. 
Circles denote the experimental data by 
Carpenter [10]. Solid curves denote the model 
predictions of Eqs. (1) and (2). 
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Figure 7. Simulated degree of cure for the 

Hercules 3501-6 resin system cured 
under constant rate of heating 
conditions indicated. 
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Figure 8. Chemovi scosi ty of the Hercules 

3501-6 resin system cured under the 
cure cycle shown. Circles denote 
the experimental data by Carpenter 
[10]. Solid curve was calculated by 
the model, Eqs. (1) and (2). 
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Figure 9. Simulated degree of cure data for the 

Hercules 3501-6 resin system cured under 
the cure cycle included. 
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